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Heavy-light physics using NRQCD-staggered actions 
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One leading source of uncertainty in the lattice computation of semi-leptonic form factors in B decay, and to 
a lesser extent B meson decay constants, comes from the extrapolation of the light quark mass to the physical 
up/down mass. This talk presents first explorations of simulating the light quark with staggered fermion actions, 
which are free of the spurious zero modes that effect Wilson-like Dirac operators and are less costly. Methods for 
fitting to Euclidean-time NRQCD-staggered meson propagators are discussed, and some preliminary spectrum 
results are presented. 



1. INTRODUCTION 

Most of the research studying heavy-light 
mesons on the lattice has focused on the formula- 
tion of heavy quarks; however systematic effects 
due to the light quark are also important to inves- 
tigate, especially in decays such as i? — > Triv. So 
far Wilson-like actions have been used exclusively 
for the light quark action, and the appearance of 
"exceptional configurations" can limit how small 
the light quark mass can be set. 

This work reports on progress applying the fa- 
miliar staggered action to the study of heavy- light 
physics, using the usual NRQCD action for the 
heavy quark. Although there are a few compli- 
cations, such as opposite parity contributions to 
correlators, this method may prove cheaper and 
more accurate than Wilson fermions since stag- 
gered fermions have only 1 spin component and 
their cutoff errors begin at O(a^) with straightfor- 
ward tree- level improvement to 0{a'^). Below we 
give a brief description of the method combining 
the two actions, followed by preliminary spectrum 
results. Observations with improved actions are 
also made. 

2. METHOD 

It is well-known that the naive discretization of 
the QCD Lagrangian 

£^ = l'(7.V + m)^', (1) 
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where 

- Ul{x-afi)-^{x-a(L)], (2) 

describes not 1 but 16 fermions due to the "dou- 
bling" symmetry. It is convenient to transform 
the fields 

■^{x) rj(a;)$(x), ■^■(x) ^ '^{x)9){x) (3) 

with Vl{x) = (--yM^^n/a to obtain a spin- 

diagonal Lagrangian 

= V + to) X 1] $ (4) 

where a^'{x) = (_i)(^o+:ri+...+x,,-i)/Q ^ jv^^^ ^-^^ 
16-fold degeneracy can be reduced to 4-fold by 
defining fields x with one spin-component: 

<^{x) = e{x)x{x) ; (5) 

usually the c-number spinor is chosen to be con- 
stant e{x) — e. K simple identity between the \1/ 
and X propagators can be derived: 

G^{x-y) = G^{x-y)^n{x)n\y). (6) 

Operators which couple to heavy-light meson 
states can be constructed from a heavy field Q{x) 
and the naive "^{x). Q{x) may be described by a 
Wilson-like discretization or a nonrelativistic for- 
mulation. Considering local bilinear operators of 
the form: 

Jrix) = Q(x)r*(x) (7) 
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with r e {1, 7'^,7'^7'',7'^7^,7^}, we can see that 
correlators of Jr can be computed using the prop- 
agator G^: 

C{p,t) = (y^e-'P-^ Tr,[GJx,t;6,0) 

\ X 

X Tr, [ri](x,i) TGq (0,0; }^ (8) 

A forthcoming paper will show that heavy-light 
mesons in this formulation do not have the com- 
plicated flavor structure that light-light mesons 
do. Mixing between pion operators is due to 
the flavor-changing exchange of high momentum 
giuons; since the heavy quark action does not 
have the doubling symmetry and describes only 
one flavor, absorption of a gluon with momentum 
component pi = ix j a would increase the heavy 
quark energy by 7r^/(2Afa^), a large factor. As 
a result, flavor-changing contributions to heavy- 
light Greens functions are suppressed. 

3. PRELIMINARY RESULTS 

An exploratory calculation was performed on 
an 8"^ X 20 lattice. 200 gauge configurations were 
generated using the tree-level, tadpole-improved 
Liischer-Weisz action within the quenched ap- 
proximation. The coupling B/^q = 1.719 cor- 
responds to 1/a = 0.8 GeV. The heavy quark 
propagators were evolved using the NRQCD ac- 
tion of Ref. 01 with bare mass oMq = 6.5 which 
roughly corresponds to the bottom quark mass. 
The light quark propagators were computed with 
several bare masses, but we focus on amo = 0.18 
which roughly corresponds to the strange sector. 
We use only local sources and sinks in this work. 

Since the NRQCD propagator requires only an 
initial condition, the heavy-light correlators are 
not periodic in t, and amplitudes and energies 
can only be extracted from (^(p, t) by fitting the 
data with t S [0,9]. Since the opposite parity 
state contributes to the correlator with a factor of 
(— 1)*+^, we must include at least two terms in the 
fit. In fact, a third is necessary since it appears 
that the first excited state is lower in energy than 
the parity partner. Thus it is desirable to utilize 
the constrained curve fitting discussed in Ref. B] . 



We extract energies and amplitudes by fitting 
the data to the form 

G(p,i) =^oe~^"* + (-l)*+iAie-^i* 

+ ^ (_l)n(t+l)^^^g-(AiJ„+AB„_. + ...)t (g-) 

n=2 

where Ai^n = _E„ — En-2- Gaussian priors were 
used for the fit parameters, with widths as fol- 
lows: SAq,i = 50%, 6Eo^i = 25%, M„>2 = 100%, 
and d{AEn>2) = 75%. The prior means for the 
ground state energies and amplitudes were se- 
lected based on effective mass plots, excited state 
amplitude means were equal to or less than the 
ground state amplitudes, and a rough guess was 
made for the excited state splittings with the as- 
sumption that AEn is roughly independent of n. 
We use all timeslices from to 9 in the fit. Most of 
our fits were done with = 3, however fits with 
A = 4 — 8 gave similar results but with increased 
uncertainties (for the priors given). When possi- 
ble, we also performed a standard unconstrained 
fit and found good agreement. 

Preliminary results are summarized in Fig. ^. 
The IS" and IP masses agree with state-of-the- 
art NRQCD/improved- Wilson calculations ||,^ 
which used finer lattice spacings. Note that the 
J = and 1 P-wave masses are automatically 
extracted as the parity partner states {Ei in (||)) 
of the Bs and B*, respectively. The 25" masses 
we extract from the AE2 in (^) are lower than 
from Ref. [|. 

4. IMPROVED ACTIONS 

Significant progress has been made in the past 
few years in improving the staggered fermion ac- 
tion, suppressing the large O(a^) flavor-changing 
interactions mediated by large momentum giu- 
ons We studied several staggered actions: 
the original 1-link action, the Naik action ||], the 
full O(a^) improved action [Q, and the "flavor- 
improved" action which is the full action minus 
the Naik term. 

The improvement of the pion dispersion rela- 
tion can be measured through the speed-of-light 
parameter c(|p|^): 



'■{m = {E'{p)~E^Q))M 



(10) 
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Figure 1. Mass splittings in the B spectrum, 
from either the 0~ or the spin averaged IS states. 
Squares mark our preHminary results, and dia- 
monds are taken from Ref. The horizon- 
tal lines correspond to experimental observations 
of the B* — Bs splitting and the spin-averaged 
B,(1P) - BsilS) splitting. 



For ap = (27r/8)fc (all permutations of k) 
we found the following preliminary values for 



action 


fc= (0,0,1) 


k = (0,1,1) 


1-link 


0.651(14) 


0.629(11) 


Naik 


0.882(14) 


0.908(36) 


flavor 


0.721(48) 


0.710(37) 


0{a') 


0.894(26) 


0.924(24) 



The Naik term is clearly responsible for improv- 
ing the pion dispersion relation. 

A good fit to the heavy-light correlators could 
not be obtained with the Naik and O(a^) actions. 
We believe this is due to the unphysical negative 
norm states present as a consequence of the 3- 
link coupling in the t direction. This hypothe- 
sis is supported by two further observations: the 
of the fits decreased when the priors allowed 
for some terms to have negative amplitudes, and 
good fits were obtained using a fifth action which 
was fully O(a^) improved in the spatial directions, 
but had no Naik term in the t direction. It is 
an open question as to why the temporal Naik 



term affects heavy-light correlators much more 
than light-light correlators. 

5. OUTLOOK 

We have demonstrated the feasibility of using 
staggered fermions as the light quark action in a 
heavy-light system. The Bayesian method of con- 
strained curve fitting is a useful, and in a few in- 
stances, necessary tool for fitting correlators on 10 
time-slices to 3 or more exponentials. The spec- 
trum on a coarse lattice agrees for the most part 
with other existing results, although the Bs{2S) 
mass extracted from our fits does not. Calcula- 
tion of the non-strange B spectrum is in progress. 
Once the appropriate renormalization constants 
are computed, this method will be useful for com- 
puting decay constants and semi-leptonic form 
factors. 
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